Anomalous quantum chaotic behavior in nanoelectromechanical structures 
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It is predicted that for sufficiently strong electron-phonon coupling an anomalous quantum chaotic 
behavior develops in certain types of suspended electro-mechanical nanostructures, here represented 
by a thin cylindrical quantum dot (billiard) on a suspended rectangular dielectric plate. The defor- 
mation potential and piezoelectric interactions are considered. As a result of the electron-phonon 
coupling between the two systems the spectral statistics of the electro-mechanic eigenenergies ex- 
hibit an anomalous behavior. If the center of the quantum dot is located at one of the symmetry 
axes of the rectangular plate, the energy level distributions correspond to the Gaussian Orthogonal 
Ensemble (GOE), otherwise they belong to the Gaussian Unitary Ensemble (GUE), even though 
the system is time- reversal invariant. 

PACS numbers: 85.85.+j, 05.45.Mt, 73.21.-b 



The possibility of engineering devices at the nano and 
micro scales has openned a great avenue for testing fun- 
damental aspects of quantum theory, otherwise difficult 
to probe in natural atomic size systems. In particu- 
lar, mesoscopic structures have played an important role 
in the experimental study of quantum chaos Q] , mainly 
through the investigation of the transport properties of 
quantum dots Q and quantum well structures 0] in the 
presence of magnetic field. However, some hard to con- 
trol characteristics of such structures can prevent the full 
observation of quantum chaotic behavior. For instance, 
the incoherent influence of the bulk on the electronic dy- 
namics hinders the observation of the so called eigen- 
state scars 4] in quantum corrals Q. Also, Random Ma- 
trix Theory (RMT) predictions to the Coulomb blockade 
peaks in quantum dots may fail due to coupling to the en- 
vironment Alternatively, suspended nanostructures, 
due to their particular architecture, are ideal candidates 
for investigating as well as implementing coherent phe- 
nomena in semiconductor devices |7|, |8j- In special, to 
understand in a controlled way how phonons influence 
electronic states and affect the system dynamics, possi- 
bly leading to a chaotic behavior. Such point is of prac- 
tical relevance since it bears the question of stability of 
quantum computers [UEi> whose actual implementation 
could be prevented by the emergence of chaos [Tl| . 

A remarkable characteristic of the quantum chaotic 
systems is the manifestation of universal statistical fea- 
tures that occur irrespective to their physical nature (e.g., 
the energy spectra of spinless particles). According to 
RMT the resulting spectra for systems with and without 
time-reversal invariance (TRI) are typically described by 
random matrices of the Gaussian Orthogonal Ensemble 
(GOE) and Gaussian Unitary Ensemble (GUE), respec- 
tively. This property was conjectured by Bohigas et al. 
[l2| and has been firmly established by theoretical and 
experimental examination [Hll3||. However, there are ex- 
ceptions to this rule, which consist of the special class 
of time-reversal invariant systems with point group irre- 
ducible representations that can exhibit the GUE statis- 



tics |3 01- The family of systems presently known 
to show this phenomenon is formed by billiards having 
threefold symmetry, which have been experimentally im- 
plemented [Tfjl IT?) in classical microwave cavities. 
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FIG. 1: a) Schematic electro-mechanical nanostruture: thin 
circular quantum dot on the surface of a suspended dielec- 
tric plate, b) Four different positions for the center of the 
quantum dot on the plate surface. 

In this letter we predict that certain electro-mechanical 
nanostructures can also exhibit the GOE as well as 
the GUE spectral statistics both under TRI condi- 
tions. We consider a thin cylindrical quantum dot (bil- 
liard) suspended on a rectangular dielectric nanostruc- 
ture (phonon cavity), as depicted in Fig. la. As the 
electron-phonon interaction is introduced between the 
electronic states, restricted to the circular quantum dot, 
and the phonon modes of the rectangular cavity an in- 
teresting interplay of the system's two relevant spatial 
symmetries takes place. In this respect, although hav- 
ing different dynamics, the present system resembles the 
much studied Sinai billiard (see, e.g., [J]), since in the 
later case chaos arises from the mismatch between the 
same two symmetries. It will be shown that for suffi- 
ciently strong electron-phonon coupling, the energy level 
distribution of this nanostructure exhibits the GOE spec- 
tral statistics if the center of the quantum dot lies in one 
of the symmetry axes of the phonon cavity, whereas the 
GUE statics occurs if center of the dot is located outside 
the symmetry axes (refer to Fig. lb). We also investigate 
the influence of material and geometrical parameters on 
the unfolding of chaos, indicating the conditions for its 
experimental observation. 



2 



The boundaries of the suspended dielectric plate can 
be either clamped (C) or free (F), corresponding to the 
Dirichlet or Neumann conditions, respectively. To ob- 
tain the vibrational modes of the plate we use the Clas- 
sical Plate Theory (CPT) approximation [l^, which is 
well suited for quasi two-dimensional mechanical sys- 
tems. The vibrations of the cavity are thus described 
by a vector field whose cartesian components are 



9W „. . dW 

U{x,y,z,t) = -z— , V{x,y,z,t) = -z— 

W(x,y,t) = J2 A m n X m (x)Y n (y)e- iut 



(1) 



In Eq. W(x, y, t) is written in terms of the one- 

dimensional transverse modes X m (x) = sin[k m x] + 
sinh[fc m a;] + £ m {cos[fc m a;] + cosh[/c m a;]}, and likewise for 
Y n . The modes X m and Yn are the solutions of the 
Bernoulli-Euler equation [l^, [Hi under the appropriate 
boundary conditions. The kinetic and strain energies for 
the cavity are 
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where p 2D , v and D represent the two-dimensional den- 
sity, the Poisson constant and the rigidity constant of 
the material, respectively. The Rayleigh-Ritz method 
is then used to obtain the coefficients Ay of Eq. QJ, 
by applying the condition dU/dAij = to the energy 
functional U = f dxdy[K.(x,y) — V(x,y)] 0|. Having 
obtained the vibrational eigenfrequencies (uj a ) and the 
eigenmodes (A^) of the dielectric cavity, an arbitrary 
deflection field is written as (Q a (t) = Q a exp[— iu>t\) 

u = J2lQ a (t)+Q* a (t)] [U a (r) i + V a (r) j + W a (r) k 

a 

(4) 

To provide the same level of description to the elas- 
tic and electronic degrees of freedom of the electro- 
mechanical nanostructure, we quantize the deflection 
field of Eq. So, we define the operator (a) that 

creates (annihilates) a mechanical mode as 

4<<> = ^QU0->/=±=>>. (5) 

Thus the quantum vibrational field that interacts with 
the electrons in the circular quantum dot is given by 



M*) + 4(*)] 

^2 Vp Lu a /h 



U a (r)i+V a (v)3 + W a (v)k 



The electrons, in the free electron gas approximation, 
occupy the states (k = (±Z, u); I = 0, 1, 2, . . .) 



J; (ai,^) exp[±iW] 
y/TrR\Ji +1 (ai u )\ 



/2 . /7TZ\ 



Here, ai v is the v-th root of the Bessel function of order 
I, R is the radius and d the width of the quantum dot. 

At low temperatures only the long wavelength acous- 
tic modes are excited and the phonon cavity can be 
treated as a continuum elastic medium. Then, the rel- 
evant electron-phonon interactions are due to the defor- 
mation (DP) and piezoelectric (PZ) potentials. The DP 
coupling operator is written as C_opA(r), where Cdp is 
the deformation potential constant and A(r) = Vu(r) is 
the relative volume variation. The DP electron-phonon 
Hamiltonian is, therefore, 



H 



DP 



h v DP 
Ci>rd— E ^bl[ai + a a ]b K ,(8) 



2Vp j 

V$L = I 4>l^{U a ,V a ,W a )ci> K dv 
Jv 

= - J2 A ^n I z [X'; n Y n + X m Y^\ 4> K dr. 

In Eq. ©, b' K (b K ) is the electronic creation (annihilation) 
operator and the integration is performed on the domain 
T>, defined by the boundary of the quantum dot billiard. 

For a piezoelectric (PZ) nanostructure of cubic crys- 
tal symmetry, the electric field produced by the cavity 
vibrations is E = [(— 2gx4,/e) e xy ] k [2£j, where gu and 
e xy are elements of the piezoelectric and strain tensors, 
respectively, and e is the dielectric constant. In the CPT 
approximation e xz = e yz = 0. Thus, the piezoelectric 
potential is (A = d (2z - d) /2) 



2Vp-J-J aa + < ] ^J- (9) 



The PZ electron-phonon Hamiltonian is, therefore, 



rPZ 



V K f ,a.K v 

= Y, A ™n I A(z) 4>* R .X' m Y^K dr. 
Jv 



(10) 



(6) 



The electro-mechanic eigenstates are thus obtained 
from the diagonalization of Hamiltonian H = Hq + 

H DP + H PZ written in the basis tyg- |0)} 

of eigenstates of H = J2 K E K b\.b R + J2 a ("a + |)^ a . 
N is the number of phonon modes included in each basis 
state. It can be verified that H is time reversal invariant. 

We investigate the dynamical behavior of the electro- 
mechanical nanostructure by performing a statistical 
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analysis of its spectrum, in terms of: (i) the eigenenergies 
nearest neighbor spacing distribution P(s), measured in 
units of mean spacing energy; and (ii) the average spec- 
tral rigidity A3 (J), which provides information about the 
correlation between energy levels within a normalized en- 
ergy interval of length I (for technical details see, e.g., 
Refs. For all the analysis we use the first 2500 

states, out of 15000 basis states, which suffice to produce 
very good spectral statistics. 

From exhaustive calculations we found that for suf- 
ficiently strong electron-phonon coupling the observed 
chaotic behavior of the nanostructure proved to be quite 
robust with respect to variations of geometrical param- 
eters, boundary conditions and basis size. Thus, for 
the sake of clearness, throughout this work we con- 
sider the representative case of a quantum dot of ra- 
dius ii = 450 nm and thickness d — 5/5 on the sur- 
face of a square phonon cavity of sides L — 1 am and 
width S = 40 nm. We also assume the clamped (C) 
boundary conditions for the four edges of the cavity. It 
has been verified that the results are similar for the dif- 
ferent boundary conditions, however the clamped only 
(X1X2Y1Y2 ={CCCC}) plate demonstrates more clearly 
the roles of the deformation potential and piezoelectric 
interactions. Nonetheless, the variation of some material 
parameters affects the chaotic behavior, for instance, the 
chaos is stronger for softer phonon cavities and for larger 
electronic effective masses. In the calculations we use the 
material parameters corresponding to an AlAs dielectric 
phonon cavity and an Alo.5Gao.5As quantum dot. This 
choice takes advantage of the very small lattice param- 
eter mismatch of the interface as well as the large effec- 
tive mass of the X valley of AlGaAs |2l|. In addition, 
to set the necessary strength of the electron-phonon cou- 
pling we multiply the corresponding interaction poten- 
tials by an arbitrary factor (3. Finally, we notice that the 
observed chaotic behavior does not result from the cou- 
pling a regular system (the quantum dot) with an already 
chaotic system (the phonon cavity). Indeed, by analyz- 
ing the spectral statistics of the vibrating plate alone we 
find a rather regular behavior. 

The interplay between the cylindrical and rectangular 
symmetries, through the electron-phonon coupling, de- 
stroys the all the geometrical invariances, except for the 
reflection symmetries if the center of the quantum dot is 
located on a symmetry axis of the plate. This scenario 
corresponds to the points A, B and D of Fig. lb. On the 
other hand, for the center of the quantum dot located at 
C, no symmetries remain. A slight displacement of the 
quantum dot is enough to generate a chaotic behavior, 
thus the relative coordinates used in the calculations are: 
A = (0,0), B = (0.05,0.05), C = (0.05,0.025), and D = 
(0.05,0). Figure |21 shows, for the DP interaction only and 
f3 = 10, P(s) and A^(l) for these four cases. In the case 
A, the spectral statistics indicates a somewhat regular 
dynamics |22( , but in B the occurrence of quantum chaos 
is clear and the level distributions are well described by 
the predictions of GOE random matrices. The same oc- 



curring for the equivalent case D. The more interesting 
result, however, is obtained for C, for which the statis- 
tics belongs to the GUE class, although the system is 
time-reversal invariant. 




FIG. 2: The level statistics for the nanostructure with the 
DP interaction only and the quantum dot centered at posi- 
tions A, B, C and D, Fig. lb. The left (right) panel represents 
by histograms (+ symbols) the numerically calculated P(s) 
(Ag(i)) distribution. The curves indicate the expected be- 
havior for regular (solid), chaotic GOE- type (dashed), and 
chaotic GUE-type (dot-dashed) systems. Here /3 = 10. 

This effect was first predicted by Leyvraz et al. 0] 
and observed experimentally pH Il7j in microwave bil- 
liards with only the threefold symmetry. In such case 
there are two classes of eigenstates: real singlets and 
complex conjugate doublets, which present the GOE and 
GUE statistics, respectively. In our case the electron 
states, Eq. J7J, naturally provide the necessary complex 
representation through the angular momentum quantum 
number /. Note that I = and I — ±1,±2, ... play, 
respectively, the role of singlet and doublet states. The 
coupling between electron and phonon systems is respon- 
sible for generating chaos and, according to Hop and 
Hpz, the electron interacts with phonon modes a, each 
of definite parity. If the dot is centered at a symmetrical 
locus, e.g., B or D, the parity makes the Hamiltonian ma- 
trix real (up to a global phase) and the GOE statistics is 
obtained. When the reflection symmetry is broken, as in 
C, the electronic angular momenta are coupled and the 
Hamiltonian is complex, exhibiting the GUE statistics (a 
detailed analysis will appear elsewhere j2^|). The same 
effect is obtained when the boundary conditions of the 
cavity are changed, which can generate, for instance, the 
GUE statistics in locus D for the {FFFC} nanostructure. 

The dependence of the spectral rigidity A 3 (/) on the 
electron-phonon coupling strength is illustrated in Fig. 
For the situation C and again considering only the DP 
interaction, we take (3 = 1, 3, 5, and 10. As f3 increases, 
the calculated statistics gradually converge to the GUE 
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prediction. Notice that the numerical results are never 
well fitted by the GOE case. The inclusion of more basis 
states does not change the observed results significantly, 
however, the effect is favored by changing some material 
parameters, for instance, by using softer phonon cavities, 
thiner plates, and including the PZ interaction. 




FIG. 3: The calculated spectral rigidity for case C and 
equal to: 1 (open circle), 3 (filled circle), 5 (x) and 10 (+). 
Once more the curves represent regular (solid) , GOE (dashed) 
and GUE (dot-dashed) type of systems. 

Finally, examination of equations © and H10|) shows 
that, individually, both the DP and PZ interactions pre- 
serve the reflection symmetry of the matrix elements with 
respect to the {x,y} axes and the diagonal axes. However, 
when acting together the reflection invariance is broken 
and the spectrum statistics of loci B and D must change 
from GOE to GUE. Fig. 0] confirms this effect by show- 
ing the P(s) distribution for the quantum dot at D, with 
the DP as well as PZ interactions included with (3 = 10. 
The agreement with the GUE statistics is excellent, in 
contrast to case D of Fig. [3 Because the AlGaAs alloy 
is a weak piezoelectric material, the DP coupling shows a 
stronger effect in promoting the chaos, whereas the main 



action of the PZ interaction (in the presence of DP) is to 
break the system's spacial symmetry. 

To conclude we briefly describe results obtained for 
other types of suspended structures , w hich will be in- 
cluded in a forthcoming publication 23]. Regarding the 
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FIG. 4: The nearest neighbor spacing distribution for case D, 
where the DP and PZ interactions are included and f3 = 10 
for both. The curves correspond to GOE (dashed) and GUE 
(dot-dashed) statistics. 

phonon cavities, we investigated different boundary con- 
ditions, in particular {FFFF} and {CFCF}. The con- 
clusions are essentially the same and the ideas presented 
here can be extrapolated to these distinct cases. We also 
considered a quantum dot of rectangular symmetry. In 
such case chaotic behavior was observed only for very 
particular asymmetric configurations, but never resulting 
in GUE statistics. This aspect shows the importance of 
the interplay between the circular symmetry of the elec- 
tronic states with the rectangular symmetry of the cavity 
phonon modes and, therefore, the fundamental relevance 
of the architecture of the nanostructures. 
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